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Abstract 



We study the effect of SU (3) flavor breaking on high density quark matter. 
We discuss, in particular, the effect a non-zero electron chemical potential 
and a finite strange quark mass. We argue that these perturbations trigger 
pion or kaon condensation. The critical chemical potential behaves as He ~ 
^mrUsA/pp and the critical strange quark mass as rus ~ m^^^A'^^^, where 
m is the light quark mass, A is the gap, and pp is the Fermi momentum. We 
note that parametrically, both the critical /ig and m1/{2pp) are much smaller 
than the gap. 
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I. INTRODUCTION 



The study of hadronic matter in the regime of high baryon density and small temperature 
has revealed a rich and beautiful phase structure One phase which has attracted 

particular interest is the color-flavor locked (CFL) phase of three flavor quark matter [Q. 
This phase is expected to be the true ground state of ordinary matter at very high density 
State of the art calculations are not sufficiently accurate to predict the critical density 
of the transition to CFL matter with any certainty. Current estimates typically range from 
Pcrit ~ (3 — 6)po5 where po is the saturation density of nuclear matter. An exciting prospect 
is the possibility to put experimental constraints on the critical density from observations 
of neutron stars. Several proposals have been made for observables that are characteristic 
of different superfluid quark phases, and attempts are being made in order to include these 
phases in realistic neutron star structure calculations P-[TO|]. 

Initial work on the superfluid phases of QCD focussed mostly on idealized worlds with 
Nf flavors of massless fermions and no external fields. But in order to understand the matter 
at the core of real neutron stars the effects of non-zero masses and finite chemical potentials 
clearly have to be taken into account. The first study of the effects of a non-zero strange 
quark mass on CFL quark matter was carried out in |Tl|,[12|. The main observation in this 
work was that a finite strange quark mass shifts the Fermi momentum of the strange quark 
with respect to the Fermi momentum of the light quarks. If the mismatch between the 
Fermi momenta is bigger than the gap then pairing between strange and non-strange quarks 
is no longer possible. The transition from CFL matter to quark matter with separate 
pairing among light and strange quarks (2+lSC) is predicted to occur at ~ -/ppK. 
Alford et al. observed that in the vicinity of this phase transition we expect to encounter 
inhomogeneous BCS phases [jl3[ analogous to the Larkin-Ovchinnikov-Fulde-Ferell (LOFF) 
phase in condensed matter physics p^^lGj. In the LOFF phase Cooper pairs have non-zero 
total momentum and as a consequence, pairing is restricted to certain regions of the Fermi 
surface. 
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In the present work we analyze CFL matter for strange quark masses and chemical 
potentials below the unlocking transition |T^. We will argue that in this regime CFL matter 
responds to the external "stress" by forming a Bose condensate of kaons or pions |18[ . This 
effect can be understood as a chiral rotation of the CFL order parameter. Superfiuid quark 
matter composed of only two flavors is characterized by an order parameter {e'^^'^u'^C'y^d'^) 
which is a flavor singlet |ll9|-pl[] . This order parameter is "rigid" and superfluidity has to be 
destroyed in order to create a macroscopic occupation number of charged excitations |^ . 
CFL matter, on the other hand, is characterized by an order parameter which is a matrix 
in color and flavor space 

where i,j labels flavor and a,b labels color indices. We can introduce a chiral field S that 



characterizes the relative flavor orientation of the left and right handed condensates ||23|. In 
the vacuum S = 1, but under the influence of a perturbation S may rotate. Because S has 
the quantum numbers of pseudoscalar Goldstone bosons, such a rotation corresponds to a 
macroscopic occupation number of Goldstone bosons. 

There is an even simpler way to explain the phenomenon of kaon condensation in super- 
fluid quark matter, see Fig. |l]. Here we concentrate on the effect of a non-zero strange quark 
mass. A non-zero quark mass shifts the energy of strange quarks in the vicinity of the Fermi 
surface by ~ m1/{2pp). In normal quark matter this leads to the decay s ^ u + e~ + (or 
s ^ u + d + u). This decay will reduce the number of strange quarks and build up a Fermi 
sea of electrons until the electron chemical potential reaches ~ m^/(4pi?). In superfiuid 
quark matter the system can also gain energy m1/{2pF) by introducing an extra up quark 
and a strange hole. This process appears to require the breaking of a pair and therefore 
involve an energy cost which is of the order of the gap A. This is not correct, however. 
An up, down-particle/strange-hole pair has the quantum numbers of a kaon. This means 
that the energy cost is not A, but mx ^ A. The CFL vacuum can decay into or 
collective modes via processes like {ds){du) + e" + z/g or — > {'us){du). 



This paper is organized as follows. In section |T| we present general arguments for the 
existence of kaon and pion condensates in high density matter with broken flavor symmetry. 
In section |T| we strengthen these arguments by performing an explicit matching calculation. 
In section |IV| we provide a different perspective on our results by using linear response theory. 

II. THREE FLAVOR QUARK MATTER AT Ms / AND / 

In order to study QCD at high baryon density it is convenient to use an effective de- 
scription that focuses on excitations close to the Fermi surface. Two effective descriptions 
of this type are available. The first effective theory is valid for excitation energies below the 
Fermi momentum pp-, while the second one applies to excitation energies below the gap A. 
The coefficients that appear in these effective theories can be worked out using matching 
arguments. In the first stage we match the microscopic theory, QCD at finite baryon density, 
to an effective theory below pp- In the second step, we match this effective description to 
an effective theory involving Goldstone modes. 

The QCD Lagrangian in the presence of a chemical potential is given by 

£ = ^ (^i;? + /i7o - /ieQ7o) ^ - i'LMijR - i^RM^L - \g%G%, (2) 

where M is a complex quark mass matrix which transforms as M — > LMR^ under chiral 
transformations {L,R) G SU{3)l x SU{3)ji, Q is the quark charge matrix, fi is the baryon 
chemical potential and /ig is (minus) the chemical potential for electric charge. As usual, 
we treat M as a (spurion) field in order to determine the structure of mass terms in the 
effective chiral theory. Once this has been achieved, we set the mass matrix to its physical 
value M = diag{mu, ma, rris). 

The quark field ip can be decomposed a.s = ip^ + ip- where ip± = |(1 ± a ■ p)ip. The 
tp^ component of the field describes quasi-particle excitations in the vicinity of the Fermi 
surface. Integrating out the ip- field we get |l2^-|26| 
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+ [r ^ L,M ^ M\Q ^ Q^) + . . 



(3) 



where = |p| — /i, = (1, v) with = p/p, A is a parameter that controls mass corrections 
to the gap and i, j, . . . and a,b, . . . denote flavor and color indices. In order to perform 
perturbative calculations in the superconducting phase we have added a tree level gap term 
ipL^fiCAipL^Ft ill the free part of the Lagrangian and subtracted it from the interacting part 
(not explicitly shown). The magnitude of A can be determined self consistently order by 
order in perturbation theory. In the normal phase both A and A vanish. In this case, only 
the first mass term in (^) contributes. 

We observe that at 0{l/pp) flavor symmetry breaking due to a chemical potential for 
charge is indistinguishable from symmetry breaking due the quark mass matrix. Indeed, up 
to terms suppressed by additional powers of (A/pp), (p/pf) or {tti/pf) the Lagrangian 
is invariant under the time dependent flavor symmetry (from now on we drop the subscript 
"+") 

(-/ieQ^ - ^) - R{t) [-^^eQ^ - ^) RKt) + iR{t)doRKt), (4) 

where L{t) and R{t) are left and right-handed time- dependent flavor transformations. 

For excitation energies below the gap A we can use an effective theory that includes only 



the pseudo-Goldstone bosons [23,26-25]. The scale of the momentum and energy expansion 
in this theory is set by the gap A. Taking into account the symmetries discussed above we 
see that a generic term in the effective lagrangian has the form 



This equation implies that the A^'th order term in the effective lagrangian is given by the most 
general SU{3)l x SU{3)ii invariant term constructed from the chiral field S and containing 
n covariant time derivatives, m spatial derivatives, p powers of M^, and q powers of fieQ 
such that N = n + 'm + p + q. We note that mass terms are suppressed by either M'^ /p\ 
or MMy {ppA). Terms of the form M'^/p\ contain the quark mass matrix in the flavor 
anti-symmetric combination shown in the gap term in Eq. (|^). 
The leading terms of the effective Lagrangian take the form 



Here S = exp(27r°A"//^) is the flavor octet chiral field and the SU{?>)a generators are 
normalized as TrfA^'A^] = 25°"^ . We have not displayed the flavor singlet part of the effective 
lagrangian. The first term in Eq. (j^) is invariant under the approximate symmetry Eq. 
because of the presence of the covariant time derivative. The second term is not invariant 
under Eq. (^, but A ~ f^A'^/pj;, is suppressed by l/pp, in accordance with Eq. (|^). 

The term is not the most general term consistent with the symmetries. The structure 
of this term determined by the fact that it has to contain the quark mass matrix in a flavor 
anti-symmetric combination. O(M^) terms that are symmetric in flavor do not vanish, but 
they are strongly suppressed. We provide an estimate of these terms in App. ^ 

Despite the similarity between the effective theory for the Goldstone modes in the CFL 
phase and chiral perturbation theory in vacuum, there are important differences in the power 
counting. As usual, the contribution of loops is suppressed by powers of p/ f-jr- However, in 
the CFL phase /^r ~ ^ A which means that the suppression of loops with respect to 
tree level terms is much more pronounced than it is in the vacuum. 

More differences appear in the expansion in M. First of all, because of an approximate 
axial Z2 symmetry in the CFL phase there are no odd powers in M. In addition to that, the 
{MM^){M^M) terms can become comparable to the terms without breaking the chiral 
expansion. Indeed, as we shall argue below, this is likely to be the case for realistic values of 
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TUg and pf- There are two reasons why the {MM^){M^ M) term can become comparable to 
the term. First, the term proportional to {MM''){M'' M) gives a contribution to meson 
masses which is of the order m? /pf while the term contributes at order m/\/pF- These 
contributions are comparable if m ~ A, which is inside the regime of validity of the effective 
theory, m < y/ ApF- Second, in the realistic case where rris ^ m^, m„, the term quadratic 
in M is proportional to at least one light quark mass, while the term quartic in M contains 
terms proportional to m^. 

Using (P) we can easily compute the masses of the Goldstone bosons in the CFL phase. 
At large density Lorentz invariance is broken and we identify the mass with the energy of a 
p = mode. For /ig = the masses of the flavored states are given by 

1/2 



ml - ml 
mjr± = =F — h 



2pi 



mK. 



T- 



1 2 

mt — m: 



2pi 



+ 
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-pr[mu + md)ms 
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-prmd[mu + ms) 



nl/2 



(8) 



2 2 

mt — mz 



m^KOKO — =F- 



2PF 



+ 



AA 
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m.u{md + ms) 
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The splitting between particles and anti-particles can be understood by observing that 
the crossed terms in the kinetic term of Eq. act as an effective chemical potential for 
strangeness/isospin even if /Xg = 0. We observe that the pion masses are not strongly affected 
but the mass of the and is substantially lowered while the K~ and are pushed 
up. As a result the and meson become massless if ~ m^^A^/^. For larger values 
of the kaon modes are unstable, signaling the formation of a kaon condensate. 

Once kaon condensation occurs the ground state is reorganized. For simplicity, we con- 
sider the case of exact isospin symmetry = = m. The most general ansatz for a kaon 
condensed ground state is given by 



E = 



exp {ia [cos(6'i)A4 + sin(6'i) cos(^^2)A5 

-|- sin(6'i) sin(6'2) cos(0)A6 + sin(6'i) sin(^; 



'2) sin 



mi]) 



(9) 



With this ansatz the vacuum energy is given by 



,2 ^2\ 2 



V{a) = -fl (^^^^^j sin(a)2 + (m?,)2(cos(a) - 1) j , (10) 

where (m^^-)^ = {'iA/ fl)mu,d{rnu,d + '^s) is the O(M^) kaon mass in the hmit of exact 
isospin symmetry. Minimizing the vacuum energy we obtain a = if m1/{2pp) < m^^ and 
cos(a) = {iTi^y / filjj- with ^eff = fnl/C^Pp) if /^e// > We observe that the vacuum 

energy is independent of 61,62, 4' even if a 7^ 0. This imphes that the effective potential in 
the kaon condensed phase has three flat directions. The hypercharge density is given by 



m 



ny = 1 - 71^ ' (11) 

where /Ze// = ml/C^PF)- This result is typical of a weakly coupled Bose gas [pOH SlH. We 
also note that within the range of validity of the effective theory, fi^ff < A, the hypercharge 
density satisfies ny < Ap|,/(27r). The upper bound on the hypercharge density in the 
condensate is equal to the particle density contained within a strip of width A around the 
Fermi surface. 

The symmetry breaking pattern is SU{2)i x f/(l)y — > U{1) where / is isospin and Y 
is hypercharge. It is amusing to note that this is the symmetry breaking pattern of the 
standard model. Kaon condensation is analogous to electroweak symmetry breaking with 
a composite Higgs field p2|,|5B|. We can discuss kaon condensation in terms of an effective 



field theory which only involves a complex kaon doublet $ = {K~^, K^) 

£ = [ (9o + tfieff) $^][ (5o - ^/ie//) - «)' ($^$) - A (12) 



If fieff > f^K the kaon field acquires a non-zero vacuum expectation value ($) = (0,f) and 
the SU{2) X [/(I) symmetry is broken to U{1). From ([T2|) we get v = {j^lff — ijmP^Y) / {2\) . 
We can fix A by comparing the amplitude of the kaon field to the result obtained from the 
chiral theory. We find A = {mPj^Y /{2fl). 

In weak coupling the coefficients of the effective Lagrangian can be computed and more 
quantitative statements about the onset of kaon condensation can be made. The gap is 



given by Q-g^J 



A = 5127r^2"i/3^2/3)-5/V^"' exp (^^^ . (13) 

The pion decay constant has been computed to leading order in (a factor 2 

discrepancy in the hterature will be resolved in section [ITl|) 



21-81og2 

18 27r2- ^ ^ 



There is also disagreement about the value of the constant A p7|j26| , |28|j38| , |39[] . The results 



given in and ||2^ are, respectively 

^log(/./A), A=- 



A=^log{i^/A), A=-^. (15) 



Using the first of these two results a condensate forms if 

^ ( 21-81og2 j ""^^^ log W^)- (16) 
In Fig. we show the dependence of the kaon mass on for pp = 500 MeV and with A 
and calculated to leading order in perturbation theory. We observe that the becomes 
massless for — 60 MeV. There is obviously some uncertainty associated with the use 
of first order perturbation theory. An estimate of this uncertainty is provided by the scale 
dependence of the result. We have calculated rriK with g evaluated at the scale A = pp. 
Varying A between Pi?/2 and 2pp gives critical strange quark masses between 39 MeV and 
67 MeV. 

If charge neutrality is enforced we have to add the contribution of electrons to the ther- 
modynamic potential, n(E, /Xe) = flcsi^, f^e) — /^e/ll^Tr^). The ground state is determined 
by minimizing Q with respect to S subject to the condition that dQ/{dfie) = 0. In the 
isospin symmetric limit these conditions are satisfied by pure K° condensation with a as 
determined above and sin(6'i) = sin(^2) = 1- This conclusion remains valid in the case 
171^ > rriu because the light quark mass difference also disfavors condensation compared 
to condensation. 

The effect of a small electron chemical potential can also be read off from Eq. (|^). A 
positive electron chemical potential lowers the energy of negatively charged Goldstone modes 
and increases the energy of positively charged modes. 



E^± = ±fie + m^±, Ek± = ±/ie + rriK^. (17) 

A meson condensate will form when /Xg equals the mass of the lightest negatively charged 
state. Let us again consider the limit of exact isospin symmetry, m„ = = m. The 
mass of the K~ is m^- = {2^/~A/ fT,)^ymms + m1/{2pF) and the mass of the tt" is mj,- = 
V^(2^/A//7r)y^mmJ. For very small the lightest negatively charged particle is the K~ , 
but for ml/{2pp) > {V2 — l)(2^/]4//^)y^mm7 the lightest negative state is the it~. For 
negative electron chemical potentials a condensate is always favored. We should note 
that the masses of charged Goldstone bosons are modified by electromagnetic effects. The 



electromagnetic self energy in the CFL phase was estimated to be ~ OemA P0| , ^ 
At sufficiently large baryon density this effect will dominate over the O(M^) contribution 
to the Goldstone boson masses. 

III. MATCHING CALCULATION FOR THE 0{M^) TERMS 



In the weak coupling regime the coefficients appearing in the Lagrangian Eq. can be 
computed by matching to perturbative QCD. In this section we will perform the matching 
calculation for the terms in Eq. (|^). Our goal is twofold: to strengthen and illustrate 
the symmetry arguments presented in the previous section and to clarify the calculations of 

in the literature^. 

We begin by calculating the one-loop polarization functions for the zeroth component 
of left-handed flavor currents ji, right-handed flavor currents Jr and (transposed) color 
currents jj. In the limit = 0, A; — > we find 



i -i 

"22 



1 _1 

2 2 



ml, (18) 



^We thank D. Kaplan for suggesting this calculation to us. 
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where the indices A,B correspond to {jL,jR,j'^) and we have introduced the quantity 



m 



21 -81og(2) / /i^ 



D 



(19) 



18 \27r^ J ' 

which is, up to a factor g"^, the Debye mass p7| , ^ . The LL and RR components of (|T8|) 
receive contributions both from diagrams with normal propagators and from diagrams with 
anomalous propagators, see Fig. ^ The LC and RC components only receive contributions 



from diagrams with anomalous propagators [^. The overall coefficient is nevertheless ex- 
actly the same. The CC entry is twice bigger than the LL and RR entries because it receives 
contributions from both left and right handed fermions. 

The matrix (0) is not diagonal, so there is mixing between gluons and left or right handed 
flavor currents. Also, there is no mixing between left and right handed flavor currents, 
contrary to what we would expect for a system with broken chiral symmetry. These defects 
can be cured by re-summing bubble chains with intermediate gluons. In practice we only 
have to compute the two-loop contribution because higher order diagrams simply correspond 
to replacing the free gluon propagator 1/ (cj^ — A;^) with the dressed propagator 1/{uj'^ — k"^ — 
g'^m'jj). The two-loop contributions to the polarization function are superficially suppressed 
by a factor g"^ but in the limit cj, /c — > the factor g"^ in the numerator is canceled by the 
screening mass g'^mj^ in the denominator. 

Summing all bubble chains we get 



4 4 " 



A 1 

4 4 ^ 



V 



m 



D- 



(20) 



1 

We observe that fiavor and color currents are decoupled and that the mixing matrix be- 
tween left and right handed current has the form expected for a system with broken chiral 
symmetry. To leading order in g'^ there are no additional contributions to the polarization 
function in the soft limit. We can now match the result (|20|) against the low energy theory 



(21) 



11 



where the covariant derivative VqS = doT, + iWiX' — iTWr determines the couphng to left 
and right handed gauge fields Wl^r. Matching the gauge field mass terms against (^) gives 
fl = mf), which is the result of Son and Stephanov p7| , P3| , p |. 

This result can also be obtained in a different way. Since the gluon field acquires a large 
mass of order gfi ^ A it does not appear in the low energy effective theory and we should 
be able to integrate it out The matrix in (|18D has eigenvalues A = —1/2, —3/2,0 and 
eigenvectors (1,— 1)/-\/2, (1, 1, — 2)/-\/6 and (1, l,l)/v^. The vanishing eigenvalue corre- 
sponds to the generators of the unbroken SU{3)l+r+c- The one-loop polarization function 
can be matched against the following mass term for the gauge fields 



D 



We 



-{Wl + Wr-2A, 



If) 



(22) 



The gauge field mass term still has the structure \{-y'){W1 + W"^ + mixing) apparent in 



). Integrating out the gluon field Aq eliminates the second term in (p2|) and we are left 



with 



m 



D 



4 2 



{Wl - Wr)' 



(23) 



which has the structure expected from the low energy effective theory (pT]) . Matching (^3]) 
against (^I]) gives = m^, as before. The important point is that in both approaches, 
summing bubble chains or integrating out the gluon field at tree level, the mixing between 
fiavor and color currents cuts down the coefficient of the quadratic terms Wl and Wr by a 
factor of 2 and introduces mixing between left and right handed currents. 

We are now in a position to perform the matching calculation for the term in Eq. (^. 
In App. ^ we present an alternative argument based on integrating out the gauge field. We 
are concerned with a possible mass term of the form 



f 

C = -— Tr 
4 

P 

= -— Tr 
2 



(24) 
(25) 



We will determine / by computing the shift in the ground state energy proportional to 
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Tr [MM^MtM] and Tr[(MMt)2] in both QCD and in the effective theory. In the effective 
theory the shift is given by 

P 



AS 



-Tr 



{MM^){M^M) - {MM^f 



(26) 



We note that the two terms in Eq. can be distinguished even in the phase E = 1 by the 
the relative position of M and M"'". We also note that other 0(M'^) terms allowed by the 
symmetries of QCD give structures that are different from the ones that appear in Eq. (p6|). 

In the microscopic theory the shift in the vacuum energy proportional to Tt[MM''M''M] 
and Tr[(MM"'')^] comes from the graphs in Figs, ^ja) and b). The Tt[MM''M''M] term is 
given by 



AS 



m 



m 



D 



1 



2 {2pf) 



D 



:Tr 



Tr 



m 



D 



Tr 



m 



D 



{MM^){M^M) 



(27) 



and the Tr[(MMt)2] term is 

AS 



m 



D 



2 (2p^)2 



Tr 



(28) 



Matching these results against Eq. (|26|) we conclude that 

which is the result we derived in section |1| from making the time derivate covariant with 
respect to time dependent flavor transformations. 



IV. LINEAR RESPONSE 

In this section we offer a different perspective on the results discussed in the previous 
sections by using linear response theory. We shall also provide a more microscopic expla- 
nation of why the two and three flavor cases behave so differently. In the three flavor case 
the system responds to a non-zero electron chemical potential by forming a condensate of 
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collective excitations. In the two flavor case, on the other hand, the response is carried only 
by the ungapped fermions. Prom an effective field theory point of view this is simply due 
to the fact that three flavor CFL quark matter has broken chiral symmetry and the low 
energy effective description contains charged collective modes whereas the two flavor theory 
has unbroken chiral symmetry and the low energy theory contains ungapped fermions and 
neutral modes. 

In order to set the stage for the discussion of superfluid quark matter we briefly review 
the response of ordinary quark matter. The grand canonical potential of non-interacting 
quarks at zero temperature is given by 



i'l — —p = 7; > 

^ 127r2 y 



ruf 



(30) 



with kf — sj 1^} — ITT'} is the Fermi momentum and and /if the chemical potential for the 
quark flavor f — u,d, s. The quark density is given by 

asi N,e, 

It is convenient to decompose the chemical potential into baryon charge, isospin, and hy- 
percharge components 

(32) 
(33) 
(34) 

We also note that /ii — V^/jLy — — A*e acts hke a chemical potential for electric charge. We 
can now study the response of the system to an external chemical potential or a change 
in the quark masses. We begin with the flavor symmetric case m„ = ma = rus = 0. The 
isospin and hypercharge susceptibilities are 







1 






= /i + 




2V3^^' 






1 




fj'd 


= // - 




2V3^^' 










fJ'S 


= - 







dni d^Vt dny d^n /^.x 
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This result has a very simple interpretation. The change in the isospin or hypercharge density 
as a function of the corresponding chemical potential is simply given by the density of states 
on the Fermi surface. The susceptibilty ( pSf ) can also be calculated in a different way, using 
the fact that x is the flavored vector current correlation function at zero momentum. We 
have 

XI = -n/(u; = 0, k^O) = ~Jd'x {j!{x)jI{0)) (36) 

with j^{x) = ip{x)'j^^ip. The correlation function (|36| ) has a vacuum piece and a density 
dependent piece. The density dependent piece is dominated by the contribution of particles 
and holes in the vicinity of the Fermi surface. We can calculate this contribution using the 



effective theory proposed in [24,25]. We get 



where ep = Ep — /i, Ep = y/jp + m?, and n(e) is the density of states. This result obviously 
agrees with Eq. (^). 

From the grand canonical potential ( pOf ) we can also determine the response of the system 
to non-zero quark masses. The derivative of the hypercharge density with respect to the 
strange quark mass is given by 



dny 



(38) 



2 „ 12? 

This result expresses the simple fact that the number of strange quarks is depleted compared 
to the number of non-strange quarks as the mass of the strange quark is increased. Again, 
we can compute this susceptibility using diagrammatic techniques. Computing a one-loop 
graph with one insertion of /iy and one insertion of m^/(2/i) we reproduce (p^). 

When we study real physical systems we are interested in the response of the system 
subject to the constraint that certain quantities are exactly conserved. In the case of neutron 
stars, for example, we are interested in the composition of quark matter subject to the 
condition that the baryon density is fixed and the net density of electric charge is zero. For 
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this purpose we consider the thermodynamic potential as a function of the quark density 
Pq = 3pB5 the up and down quark fractions x = Pu/Pq and y = Pd/ Pq, and the electron 
chemical potential pe 

uip„ X, y, pe) = F{p„ X, y) - p^Q = ^pf {x'/' + y'/' + {1 - x - yf/' 

+n~^/'p-'/'ml{l - X - yf'} + p.p. - ^) - ^,pt- (39) 

We have neglected higher order terms in the strange quark mass as well as the mass of 
the electron. In order to determine the ground state we have to make ([39| ) stationary with 
respect to x,y,pe- Minimization with respect to x and y enforces j3 equilibrium, while 
minimization with respect to pe ensures charge neutrality. We find 

2 

TTI 

which shows that there is a small non-zero pe and a corresponding suppression of strange 
quarks with respect to light quarks even at high density. 

We would now like to study how these results are modified in superfluid phases of QCD. 
We begin with a simple toy model introduced by Rajagopal and Wilczeck ||T^. The model 
contains two quark flavors, up and down, that pair in a spin singlet state which is anti- 
symmetric in both color and flavor. The pair condensate is described by the order parameter 
{e"'^u"'C'~f5d^). Here, a,b are color indices that only take on the values 1 and 2. One may 
think of this toy model as Nf = 2 QCD where the contribution of the third, unpaired, quark 
color is ignored. Alternatively, we may think of this theory as A''^ = 2 QCD. 

We can calculate the response in the superfluid in the same way we did in the normal 
phase, using the relation between the quark number susceptibilities and the 00-component 
of the polarization function. In the superfluid phase there are two contributions, coming 
from the normal and anomalous components of the quark propagator. For the quark number 
susceptibility we get 

XB = -noo(u; = 0,fc - 0) = 4Ar, / ^ { (^g -e^l^.y - (^g^f^A^) ' ^^1) 
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where the first term is the contribution from the normal quark propagator and the second 
term is the anomalous contribution. The two contributions are exactly equal and sum up to 

which is equal to the result in the normal phase. We should note that the first term alone 
only contributes half the susceptibility in the normal phase, even though the susceptibility 
is independent of the gap and the naive A ^ limit of the first graph would seem to 
correspond to the susceptibility in the normal phase. This is due to the fact that the u ^ 
and A — > limits do not commute. This phenomenon is well known from calculations of 
the screening mass in other many body systems |^ . 



The calculation of the isospin susceptibility proceeds along exactly the same lines, only 
the isospin factors of the two diagrams are different. The isospin factor of the normal 
contribution is tr[r3r3] = 2, while the isospin factor of the second term is tr[r3r2r3r2] = —2. 
The two contributions cancel exactly and the isospin susceptibility is zero. This results has 
a simple physical interpretation. The superfluid order parameter in Nf = = 2 QCD is 
a flavor singlet and the only broken symmetry is the U{1) of baryon number. As a result 
there is only one massless state, the f/(l) Goldstone boson. This state couples to the baryon 
density and leads to a non-zero baryon number susceptibility but it does not couple to 
isospin. All states that carry isospin have energies of the order of the gap, so xi remains 
zero as long as /i/ < A. 

We can also see how the calculation of the isospin susceptibility differs in the case of 
CFL quark matter. Because of the symmetries of the CFL phase there are two types of 
quasi-particles, an SU{3) octet with gap Ag = A and an SU{3) singlet with gap Ai = 2A. 
Up to degeneracy factors the two types of quasi-particles contribute equally to the quark 
number susceptibility. We find xb = ^Sfi"^/ (27r^) which is equal to the result in the normal 
phase. The calculation of the isospin susceptibility is more complicated. We get 



(27r)4 I 6 (p§ - e2 - Ai)(p2 _ - Ai) 3 {pi - - Al)(pl - - Af) 
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(43) 



_ 1 A| 1 AgAi 1 

3 (pI -4- ^DiPo -4- Ai) 3 (pI -el- Aj^ - ej - 4A?) / 

The first term comes from particle-liole diagrams witli two octet quasi-particles wliile tlie 

second term comes from diagrams with one octet and one singlet quasi-particle. There is 

no couphng of an octet field to two singlet particles. The third and fourth term are the 

corresponding contributions from particle-particle and hole-hole pairs. The four integrals in 

m give 



7 1 _ 1 _ 41og(2) ] ^ 21-81og(2) ^ . . 

6 3 3 9 I 1 4vr2 ] 18 I 27r2 i ' \27i^ 



which should be compared to xi = 3/i^/ (27r2) in the normal phase. We observe that there is 
a partial cancellation between the normal and anomalous contributions. However, because 
of the more complicated flavor structure this cancellation is not exact. The isospin density 
induced by an isospin chemical potential is reduced by a factor ~ 3.5 compared to the 
normal phase, but it does not vanish. In linear response theory we expand around the point 
/i/ = m„ = = rUs = 0. In the real world the quark masses are non-zero and there is a 
critical isospin chemical potential ^f^^ ^ below which the isospin susceptibility vanishes. 
In order to see a threshold behavior in fij we have to resum mass corrections. This is most 
efficiently accomplished using the effective chiral description developed in section || , see 
equation ([10|) . 

V. SUMMARY 

We have studied the response of three flavor quark matter to a non-zero electron 
chemical potential and a non-zero strange quark mass. We have focussed on the regime 
^ei^ll {'^Pf) < a in which the perturbation does not destroy color-flavor locking. We have 
identified a new scale /Xg, ''^^/(^Pf) ~ y/^u,d^s{^/PF) which corresponds to the onset of 
pion or kaon condensation p6|^T|. This scale is parametrically much smaller than the gap. 
If CFL quark matter exists in the core of a neutron star it is likely to be condensed. 



Both with or without a kaon condensate there are no electrons present IT^- If CFL quark 



matter is in contact with a hadronic phase that supports a large electron chemical potential 



the surface layer is likely to be K~ or tt~ condensed |]52 |. 

These results are based on an analysis of how to incorporate fie and m^/ {2pp) in the chiral 
effective theory. Both terms enter as constant flavor gauge fields, with coefficients completely 
determined by The contribution of the m^/ {2pp) term to the Goldstone boson masses is 
of higher order in the quark masses as compared to the leading order y/mms{A/pF) term. 
It can nevertheless become dominant because the 0{m) term is suppressed by powers of 



m/nis and (A/pp)- As a consequence the 0{m?) term can cancel the 0{m) term without 
leading to a breakdown of the low energy expansion. 
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APPENDIX A: MASS TERMS INDUCED BY THE COLOR SYMMETRIC 

DIQUARK CONDENSATE 



The O(M^) mass term in gives anomalously small Goldstone boson masses of the 
order uigb ~ y/nim^i'^/pF)- We already noted that mass terms not suppressed by (A/pp) 
cannot appear at O(M^). For strange mesons the O(M^) mass term also contains an ad- 



ditional suppression factor ym/m^. Here, m is the mass of the light quarks and is the 
strange quark mass. The fact that all Goldstone boson masses are proportional to the light 
quark mass is related to the fact that the CFL order parameter is totally anti-symmetric in 
flavor. This flavor structure also leads to an accidental symmetry of the effective theory at 
O(M^). If rris = but m 7^ we find an octet of exact Goldstone bosons, even though the 
unbroken flavor symmetry is only SU{2). 

There are mass terms at O(M^) that are consistent with the symmetries of the CFL 
phase that will remove the accidental symmetry and give contributions to the kaon mass 
that are proportional to ms{A/pp). These terms are induced by the color-flavor symmetric 
gap parameter 

Ag' = As [m] + ^".^f) ■ (Al) 

The symmetric gap is consistent with the symmetries of the CFL phase but disfavored by 
the interaction. In particular, one-gluon exchange is repulsive in the color-symmetric quark- 
quark channel. In perturbative QCD, a small symmetric gap is generated by mixing with 
the primary gap parameter. We find 

where A^ is the color-flavor anti-symmetric gap parameter. 

We can calculate the contribution of A5 to the Goldstone masses using the methods of 



Beane et al. |]26|. Including the effects of both A^ and A5 we find 
A 



C = -^^^ log ( — ) (Tr(MS)Tr(MS) - Tr(MSME) + h.c.) 
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- log j (Tr(MS)Tr(MS) + Tr(MSMS) + h.c.) . (A3) 

Here, Aa,s are the flavor anti-symmetric and symmetric "anti-gaps". For the purpose of 
estimating the relative size of the two mass terms we shall assume that Aa,s — Aa,s- We 
can now calculate the correction to the charged kaon mass, 

~ — - — ^mnis — + ••• , A4 

Jtt V V^A/ mu ) 

with Aa = A^/(47r^) log(pi?/A^). Using ( [A^ ) we observe that even for ms/rriu — 20 the 
correction to the kaon mass due to the color-symmetric gap is still small. 

APPENDIX B: /iQ + MMt/(2PF) TERMS FROM INTEGRATING OUT THE 

GAUGE FIELD 



Following the discussion in section |ITl| we can also derive the O(M^) terms by integrating 
out the gauge field. This discussion will also make it clear that the M^M and MM^ terms 
enter in the effective lagrangian like gauge fields, together with flavor non-singlet chemical 
potentials. 

In this section we would also like to show how, by explicitly keeping track of the orien- 
tation of the CFL order parameter, we can determine how the chiral field S enters into the 
mass terms. This is useful because at higher order the number of independent terms in the 
chiral lagrangian quickly proliferates and it becomes more difficult to identify the diagrams 
in the microscopic theory that correspond to a given term in the effective lagrangian. 

In order to match the microscopic theory against the effective theory in the vacuum 
(S = 1) phase we calculate diagrams in the microscopic theory using the Nambu-Gorkov 
propagators in the normal CFL phase. The inverse Nambu-Gorkov propagator for the ip^ 
field is given by 
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s 



with the anomalous self energy 



-1 



V 



A po + e 



(Bl) 



(Ai)^l = -{Ar)^^ = As f(5f(5? - 6'}6' 



(B2) 



The inverse Nambu Gorkov propagator is not diagonal in color and flavor. It becomes 
diagonal in the space spanned by the 9 color-flavor matrices 

^ 

JjCA^"- 



,A\ai 



V 



(B3) 



where A° = ^^2/3 and A^, (A = 1, . . . , 8) are the Gell-Mann matrices. In this basis it is 
straightforward to compute the inverse of ( pT|) . We flnd 



S 



AB 



■iAB 



( Po + e -A^ ^ 



.2 _.2_M 



V 



-A^ PQ-t, 



(B4) 



with A^ = 2A8 for A = and A^ = -sym(A)A8 for A = 1, . . . , 8. Here, sym(y4) = 1 for the 
symmetric Gell-Mann matrices A = (1,3,4,6,8) and sym(A) = —1 for the anti-symmetric 
matrices A = (2, 5, 7). 

In order keep the dependence on S we have to perform the calculation using the anoma- 
lous self energy in the rotated vacuum 



(A^)g' = As (xf - 



ab 



As 



(B5) 
(B6) 



with X G SU (3)l and Y G SU (3)/j. The Nambu-Gorkov propagator for left handed fermions 
is diagonal in a basis spanned by the color-flavor matrices 



( 1 l\AvT\ai 



V 



^(A^Xt)^ 



(B7) 
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with a similar set of matrices ("5^)"* which diagonahze the propagator for right handed 
fermions. In the basis ( [B7| ) the fermion propagator in the rotated CFL vacuum has exactly 
the same form (P^ ) that it had in the ordinary CFL vacuum ( |B2| ) . The dependence on X, Y 
comes in when we calculate diagrams with external color or flavor currents. In that case we 
have to take matrix elements of the external current between the basis states (vl) and (vr). 

We can now calculate a one-loop diagram with insertions of MM"'" and the gauge fleld 
Aq. We find 



AS 



m 



D 



2pi 



Tr 



(B8) 



In the same way, we also calculate diagrams with insertions of M^M and Q. Collecting all 
these terms we get 

2 / , ^+ , ^ \ 2" 



m 



D 



Tr 



(B9) 



Similar to the calculation of /^r it is essential here to take into account the mixing with the 
gauge field. Without the Aq field we would conclude that there is no dependence on the 
fiavor matrices X, Y. We can now integrate out the gauge field Aq. We get 



AS 



m 



D 



4 



Tr 



MtM\ 

2PF I 



MMty 



(BIO) 



where S = XY'^ . We note that after integrating out the gauge field the vacuum energy 



( BIO ) only depends on the chiral field S and not on X and Y separately. Using = tno 
we observe that (|B10| ) contains the terms required to complete the covariant derivative in 



23 



REFERENCES 



[1] K. Rajagopal and F. Wilczek, |hep-ph/001133^ . 
[2] M. Alford, |hep-ph/0l02(M7. 



[3] T. Schafer and E. Shuryak, |nucl-th/OOlO[HU 



[4] M. Alford, K. Rajagopal and F. Wilczek, Nucl. Phys. B537, 443 (1999) jhep^ 
ph/ 98044031 . 



[5] T. Schafer and F. Wilczek, Phys. Rev. Lett. 82, 3956 (1999) ||hep-ph/981147^" . 
[6] T. Schafer, Nucl. Phys. B575, 269 (2000) ||hep-ph/ 990957^ . 



[7] N. Evans, J. Hormuzdiar, S. D. Hsu and M. Schwetz, Nucl. Phys. B 581, 391 (2000) 
|hep-ph/99103T3 . 



D. Page, M. Prakash, J. M. Lattimer and A. Steiner, Phys. Rev. Lett. 85, 2048 (2000) 
^ep-ph/00050gi. 



[9] D. Blaschke, T. Klahn and D. N. Voskresensky, Astrophys. J. 533, 406 (2000) pstra 



ph/9908334| . 



[10] M. Alford, J. A. Bowers and K. Rajagopal, J. Phys. G27, 541 (2001) ||hep-ph/0009357 



[11] M. Alford, J. Berges and K. Rajagopal, Nucl. Phys. B558, 219 (1999) ||hep-ph/9903502 



[12] T. Schafer and F. Wilczek, Phys. Rev. D60, 074014 (1999) [|hep-ph/9903503 |. 



[13] M. Alford, J. Bowers and K. Rajagopal, Phys. Rev. D 63, 074016 (2001) jEi^ 
I ph/0008208| . 



[14] A. 1. Larkin and Yu. N. Ovchinnikov, Zh. Eksp. Theor. Fiz. 47, 1136 (1964); engl. 
translation: Sov. Phys. JETP 20, 762 (1965). 

[15] P. Fulde and A. Ferrell, Phys. Rev. 145, A550 (1964). 



24 



[16] A. A. Abrikosov, Fundamentals of the theory of metals, North-Holland, Amsterdam 
(1988). 



[17; 

[18 
[19 
[20 



[21 

[22 
[23 
[24 
[25 
[26 



[27 

[28 
[29 



K. Rajagopal and F. Wilczek, Phys. Rev. Lett. 86, 3492 (2001) ||hep-ph/OOl:^IIH^ . 
T. Schafer, Phys. Rev. Lett. 85, 5531 (2000) | |nucl-th/OOOrTm |. 



D. Bailin and A. Love, Phys. Rept. 107, 325 (1984). 

M. Alford, K. Rajagopal and F. Wilczek, Phys. Lett. B422, 247 (1998) jhep^ 



ph/9711395| . 



R. Rapp, T. Schafer, E. V. Shuryak and M. Velkovsky, Phys. Rev. Lett. 81, 53 (1998) 
|hep-ph/9711396 |. 



P. F. Bedaque, |hep-ph/9910247 . 



R. Casalbuoni and D. Gatto, Phys. Lett. B464, 111 (1999) ||hep-ph/990822^ . 



D. K. Hong, Phys. Lett. B 473, 118 (2000) ||hep-ph/98l25TI]|] . 
D. K. Hong, Nucl. Phys. B 582, 451 (2000) ||hep-ph/99053^ . 



S. R. Beane, P. F. Bedaque, and M. J. Savage, Phys. Lett. B483, 131 (2000) ^epj 



ph/0002209 . 



D. T. Son and M. Stephanov, Phys. Rev. D61, 074012 (2000) [|hep-ph/9910491|| , erra- 
tum: |hep-ph/ 0004095 . 



C. Manuel and M. H. Tytgat, Phys. Lett. B479, 190 (2000) ||hep-ph/0001095[ . 



D. T. Son and M. A. Stephanov, Phys. Rev. Lett. 86, 592 (2001) ||hep-ph/0005"225 



[30] J. B. Kogut, M. A. Stephanov, D. Toublan, J. J. Verbaarschot and A. Zhitnitsky, Nucl. 
Phys. B 582, 477 (2000) ||hep-ph/OOOlTn |. 



[31] J. B. Kogut and D. Toublan, Phys. Rev. D 64, 034007 (2001) | |hep-ph/0l03?n^ . 



25 



[32] E. Farhi and L. Susskind, Phys. Kept. 74, 277 (1981). 

[33] D. Kaplan and H. Georgi, Phys. Lett. B 136, 183 (1984). 

[34] T. Schafer and F. Wilczek, Phys. Rev. D60, 114033 (1999) | |hep-ph/9906^ |. 



[35] R. D. Pisarski and D. H. Rischke, Phys. Rev. D61, 074017 (2000) ||nucl-th/99lOTJ3B 



[36] D. K. Hong, V. A. Miransky, I. A. Shovkovy and L. C. Wijewardhana, Phys. Rev. D61, 
056001 (2000) [|hep-ph/9906478| . 



[37] W. E. Brown, J. T. Liu and H. Ren, Phys. Rev. D61, 114012 (2000) ||hep-ph/9908"248 



[38] M. Rho, A. Wirzba, and L Zahed, Phys. Lett. B473, 126 (2000) ||hep-ph/9910"550[ . 



[39] D. K. Hong, T. Lee, and D. Min, Phys. Lett. B477, 137 (2000) ||hep-ph/9912"53T 
[40] D. K. Hong, Phys. Rev. D 62, 091501 (2000) ||hep-ph/0006T05| . 



[41] C. Manuel and M. H. Tytgat, Phys. Lett. B 501, 200 (2001) ||hep-ph/ 00 10274 



[42] D. H. Rischke, Phys. Rev. D 62, 054017 (2000) ||nucl-th/0003'063 



[43] K. Zarembo, Phys. Rev. D 62, 054003 (2000) ||hep-ph/0002"T23| . 



[44] V. A. Miransky, I. A. Shovkovy and L. C. Wijewardhana, Phys. Rev. D 63, 056005 
(2001) ||hep-ph/ 00091731 . 

[45] A. A. Abrikosov, L. P. Gorkov and I. E. Dzyaloshinski, Methods of quantum field theory 
in statistical physics, Prentice-Hall, Englewood Cliffs, N.J. (1963). 

[46] A. B. Migdal, Sov. Phys. JETP 36, 1052 (1973). 

[47] R. F. Sawyer, Phys. Rev. Lett. 29, 382 (1972). 

[48] D. J. Scalapino, Phys. Rev. Lett. 29, 392 (1972). 

[49] D. B. Kaplan and A. E. Nelson, Phys. Lett. B175, 57 (1986). 



26 



[50] G. E. Brown, K. Kubodera, and M. Rho, Phys. Lett. B175, 57 (1987). 

[51] H. D. Politzer and M. B. Wise, Phys. Lett. B273, 156 (1991). 

[52] M. Alford, K. Rajagopal, S. Reddy and F. Wilczek, |hep-ph/0l05(]Ug. 



27 



FIGURES 




FIG. 1. Schematic picture of weak decays in normal (a) and superfluid (b) quark matter with 
three quark flavors. We assume that initially the density of all quark flavors is the same, so that 
ep^s — ^F,ud + 1^1/ {'^Pf)- Solid and open circles show particles (p) and holes {h). In (a) a strange 
particle decays into an up quark, an electron and a neutrino, leaving behind a strange hole. In the 
left panel of (b) a strange particle decays into an up quark, a down particle-hole pair, an electron 
and a neutrino. The remaining {pp){hh) configuration has the quantum numbers of a K'^. In the 
right panel we show the decay of a strange quark into a {pp){hh) configuration with the quantum 
numbers of a K^. 
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FIG. 2. Masses of and excitations in the color-flavor locked phase. We show the 

excitation energies as a function of for = 500 MeV. The gap A = 67 MeV and the pion 
decay constant /tt = 104 MeV were determined to leading order in perturbation theory. The solid 
and dashed curve show the masses of the {K^^K^) and {K~ ,K^) states. The dotted curve shows 
the kaon masses calculated from the leading order 0{mq) term. The short dashed curve shows the 
pion masses. 
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a) 




FIG. 3. Diagrams contributing to the two-point functions of two L currents (Fig. a), one L 
and one R current (Fig. b), and one L and one color current (Fig. c). The squares denote the 
anomalous fermion self energy while the triangle denotes a resummed gluon propagator. 
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a) 




FIG. 4. Fig. a) shows the diagram in the microscopic theory which is matched against the 
term in the chiral theory. Fig. b) shows the diagrams which are matched against 
the (MMt)2 term. 
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